Proposal for entangling remote micromechanical oscillators via optical measurements 



K. B0rkje,' A. Nunnenkamp,' and S. M. Girvin' 

'Departments of Physics and Applied Physics, Yale University, New Haven, Connecticut 06520, USA 

(Dated: Received January 14, 2013) 

We propose an experiment to create and verify entanglement between remote mechanical objects by use of 
an optomechanical interferometer. Two optical cavities, each coupled to a separate mechanical oscillator, are 
coherently driven such that the oscillators are laser cooled to the quantum regime. The entanglement is induced 
by optical measurement and comes about by combining the output from the two cavities to erase which-path 
information. It can be verified through measurements of degrees of second-order coherence of the optical output 
field. The experiment is feasible in the regime of weak optomechanical coupling. Realistic parameters for the 
membrane-in-the-middle geometry suggest entangled state lifetimes on the order of milliseconds. 
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Despite the tremendeous success of quantum mechanics 
at explaining the behaviour of the microscopic world, many 
people have found it uncomfortable that macroscopic objects 
should also obey the laws of quantum mechanics. This goes 
back to the founders of quantum mechanics, and is the ori- 
gin of the famous Schrodinger cat thought experiment IT]. 
However, if all time evolution is unitary according to the 
Schrodinger equation, as in Everett's relative-state interpreta- 
tion ||2l, there is nothing that forbids counterintuitive phenom- 
ena such as superpositions of macroscopically distinct states 
f3 1. It should therefore be possible to observe quantum effects 
on arbitrarily large scales if the system is adequately shielded 
from environmentally induced decoherence. Conversely, ex- 
periments might reveal the existence of unknown sources of 
decoherence |4 ,3'|, which limit quantum mechanics to small 
scales by causing an objective wavefunction collapse and cre- 
ating a quantum-classical boundary. 

Cavity optomechanics, where mechanical oscillators are 
coupled to light or microwaves IIHSI, is a promising field for 
experimental tests of quantum mechanics at large scales. The 
motion of a micromechanical object can be cooled via radia- 
tion pressure forces 191 4161 . One group has already reported to 
have reached the quantum regime by resolved sideband cool- 
ing 1 17 1, and others are expected to soon follow. An interest- 
ing future direction for experiments would be to detect quan- 
tum entanglement |18 | between a micromechanical oscillator 
and another system, such as an optical cavity field mode 1191 
or a second mechanical oscillator ll20l - [23l . 

In this Letter, we propose a new and promising route to cre- 
ate and verify entanglement between two remote mechanical 
oscillators. Our idea is based on an optomechanical interfer- 
ometer where two optical cavities, each coupled to a separate 
mechanical oscillator, are coherently diiven in parallel. The 
drive frequency is chosen such that the oscillators are cooled 
close to the motional ground state. In short, entanglement is 
achieved by combining the optical output from the two cavi- 
ties to erase which-path information. The mechanical objects 
in our setup have no direct interaction, but are projected onto 
an entangled state through optical measurements, in contrast 
to various earlier proposals f2^^TT\. A scheme related to ours 
involving laser pulses was discussed in Ref. |23|. Similar ap- 
proaches have been successfully applied to entangle remote 
atomic ensembles [i24i i25 J as well as individual trapped ions 



ll26l . We also present a new and feasible way of detecting 
the entanglement. The experiment we propose should be re- 
alizable with present day technology in the regime of weak 
optomechanical coupling. For the membrane-in-the-middle 
geometry [16|, we estimate that the entangled states can have 
decoherence times on the order of milliseconds. This could 
have relevance to quantum information and communication 
technologies |27 , 28|. 

Consider the standard optomechanical system where the 
position of a mechanical oscillator modulates the frequency 
of an optical cavity mode. We discuss a single cavity first and 
move on to a system with two cavities later. We will have 
the membrane-in-the-middle setup 1 16| in mind, but our dis- 
cussion can apply to many other realizations. The photon and 
phonon annihilation operators will be denoted by a and c, re- 
spectively. The interaction Hamiltonian is i?i,it = hgxa^a 
where g is a coupling constant. The mechanical position op- 
erator is = Xzpf (c + ) , with x^pf being the size of the zero 
point fluctuations. Input-output theory 1.29. .30J leads to the 
quantum Langevin equations 

a = - (^'^ +iujc^ a~ igxa + '^^,ain,i , (1) 

i 

c = - (^'^ + iujuj c ~ igx^pta^ a + . 

The bare mechanical and optical frequencies are ujm and utc- 
The mechanical oscillator is coupled to a thermal bath, re- 
sulting in a nonzero linewidth 7 and a fluctuating force deter- 
mined by the operator Cin- For systems where the mechanical 
quality factor wm /7 is high, the Markov approximation gives 

[cinW,4(^')] = Kt-t') and (4(0cin(i')) = nth^(i-i')- 
Here, nth ~ k-QT/hujM, where T is the bath temperature. 
The optical linewidth is k ~ J^i where Ki is the decay rate 
through decay channel i. We imagine a two-sided cavity with 
a left (L) and a right (R) input pott, and assume that the cavity 
is driven from the left. The optical input modes then take the 
form ain,L e"*'^°*ain + and ain,R = init), where ai„ 
is a constant, wd is the drive frequency, and the vacuum noise 

operators 6 obey ^.(i), - = %<S(i-i')- 

We define the detuning between the drive and the mean cavity 
frequency as A = wd — ^ 9{x)- 

The cavity mode operator can be written as a sum of a mean 
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and a fluctuating part, a{t) = e ^'^^*{a + d{t)), where |ap 
is the mean number of photons in the cavity. Here, we fo- 
cus on the situation (d^d) <^ \d\^. In that case, Eqs. can 
be linearized and solved analytically. The effective coupling 
between the optical and mechanical fluctuations is given by 
a — gx^pfd. In the regime of weak coupling |a| <C n and 
for negative detuning A, the mechanical oscillator is approx- 
imately in a thermal state, but with renormalized parameters 
compared to the case of g = 0. The frequency Cjm is shifted 
from its bare value due to the interaction with the optical field, 
often referred to as the optical spring effect |31 1. The effec- 
tive linewidth 7 = 7 + 7opt is now a sum of the bare value 
and a contribution 7opt due to the optomechanical coupling. 
The latter is positive when the detuning A is negative, lead- 
ing to line broadening. The effective mean phonon number is 
the weighted sum t^m = (7"-th + 7opt"-opt)/7. where jiopt is 
a measure of the effective temperature of the radiation pres- 
sure shot noise ll30l l32l l33l . Below, we will set the detuning 
to A = — wm, which is optimal for cooling, and assume that 
7opt ^ 7 and riopt — k^/(4wm)'^ < 1- Note however that 



our results are also of interest when riopt > 1 (see Ref. 11341 ). 

The optical output mode on the right-hand side of the cavity 
is b{t) = y^KRa(i) — ^R{t). The optomechanical interaction 
leads to sidebands in the output which are displaced from the 
drive frequency oj-d by the mechanical frequency com- The 
sidebands originate from Raman scattering off the mechanical 
oscillator, where a photon gains (loses) energy by destroying 
(creating) one phonon. We will imagine that photons at the 
drive frequency can be completely filtered away and that the 
two sidebands can be measured independently. We will refer 
to the output at frequency ojd ± wm as the blue (red) sideband, 
and denote the output mode filtered around this frequency by 
bh{r){t)- See the Supplementary Material t34J for details. The 
width of the sidebands are given by the mechanical Unewidth 
7, so we assume the filter bandwidth A to obey 7 ^ A ^ ujm- 
The ratio between the output fluxes of red and blue photons is 

noptiriM + l)/(»^opt + l)n-M- 

A filter that removes the carrier photons at and splits the 
red and blue photons into different spatial modes might pose 
a technical challenge. One reason is that the mechanical fre- 
quency is typically in the kHz-MHz range. However, there are 
also experimental setups with um in the GHz range Il35l437l . 
Another reason is that the ratio between the fluxes of blue and 
carrier photons, given by A{gx2pf / K)^nM, is very small in the 
weak coupling limit. An alternative and more feasible way 
to achieve the filtering is through heterodyne photodetection, 
where the blue and red sidebands can easily be distinguished 
in the Fourier domain. This is discussed in detail in the Sup- 
plementary Material |34|. 

We define the degrees of second-order coherence 1 
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where steady state has been assumed, r > 0, and the in- 
dices i and j denote either red (r) or blue (b). We find that 
(''') = 9h\h(^) ~ l + e^^"^. This is what one would expect 
for thermal radiation seen through a Lorentzian filter of width 



7 1*381. The photon statistics of the red and blue sidebands is 
that of thermal radiation simply because the mechanical oscil- 
lator is approximately in a thermal state. More interestingly, 
the cross-correlators become 
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when restricting the delay time such that t ^ k^^, A~^ f34l. 
The result ([3]l has a simple physical explanation. If the ef- 
fective phonon number um is less than 1, the mechanical os- 
cillator spends most of the time in the ground state |0). It 
can gain energy and reach the excited state |1) by the cre- 
ation of a red photon. However, it is bound to return to the 
ground state quickly, through the creation of a blue photon. 
This means that conditioned on the detection of a red pho- 
ton, the probability of detecting a blue photon is high, such 
that g^^ir) becomes large. The opposite is not the case. For 
TT-M < 1, once a blue photon is detected, it probably means 
that the oscillator is now in the ground state |0) and detec- 
tion of a red photon is not particularly likely. In fact, in the 

limit ?iM 0, having detected a blue photon does not change 

(2) 

the probability of detecting a red one, such that i^^n^ (t) — > 1. 
Note that for sufficiently small riM, the classical inequality 
.(2)/^M2 ^ „(2)^n^„(2)^m I.38J can be violated. This hap- 

> {bl{t)k{t)){bi{t)bi,{t)) and 



pens when \ {bl{t)bl{t 
means that the fields 6b and 4 are entangled. This is a nec 
essary requirement for detecting entanglement between me 
chanical oscillators in the setup that we discuss next. 



Cavity 1 




FIG. 1 : (Color online) Schematic view of our proposed experimental 
setup. Combining the output from the two cavities on a beam splitter 
can create entanglement between the mechanical oscillators. This 
can be verified by measuring the photon statistics of the red (dashed) 
and the blue (dotted) sidebands. The sideband filtering can also be 
achieved through heterodyne photodetection. 



We now move on to the main part of this Letter and study 
the setup presented in Fig. [T] We consider two optical cav- 
ities, 1 and 2, each coupled to a mechanical oscillator. For 
simplicity, the cavities are assumed to be identical, such that 
wc,i ~ wc,2 = i^Ci Ki = K2 = K, etc. We imagine that 
kr ^ kl, such that most photons leave the cavity through 
the mirror on the right. Furthermore, we assume that the op- 
tomechanical coupling constants are equal. We do however 
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allow for the two mechanical oscillator frequencies to differ, 
defining ujm,i = wm and wm,2 — ^^m + but we restrict the 
difference to \S\ ^ k,ujm- The cavities are driven in parallel 
by a laser optimally detuned to A = —ujm, such that cool- 
ing to the quantum limit is achieved for both oscillators. The 
right-hand side output from the two cavities are combined on 
a 50:50 beam splitter and then filtered into red and blue side- 
bands, as discussed above. The red photons are detected at 
either the photomultiplier Ap or B^, whereas the blue photons 
are detected at photomultiplier Ab or Bb. 

Let us first discuss how this setup can lead to entangle- 
ment between the two mechanical oscillators. Assume that 
the two oscillators are identical and in the state |0,0), i.e., 
zero phonons in both. The detection of a red sideband pho- 
ton means that one mechanical excitation, a phonon, has been 
created. However, the information on which cavity the photon 
came from has been erased. This means that, conditioned on 
detecting a red photon at either Aj or Bi, the wavefunction 
of the oscillators has now collapsed onto the superposition 
(|1, 0)+e*^|0, l))/\/2. This is an entangled state, even though 
it contains only one phonon [3^i40|. The phase depends on 
whether the photon was detected at A^ or B^, and other factors 
such as optical path length differences. While this simplified 
discussion provides insight on how entanglement is created, 
the thermal baths must of course be taken into account. 

It is interesting to examine the degrees of second-order co- 
herence defined in Eq. (j2]i for this setup. Now we need to 
specify not only photon "color", but also which detectors we 

(2) 

refer to. Taking 5A.b|A ^'^ example, it is instructive to 

express this as 



(4) 



defining t' = t + t. Here, the expectation value in the 
denominator is the photon flux at detector Ab with respect 
to the state pss, which is the steady state density matrix in 
the absence of measurements. In the state p^s there is ob- 
viously no entanglement between the mechanical oscillators. 
On the other hand, the expectation value in the numerator is 
defined by (0(i'))A, = Tr(Op(i')), where, formally, p{t') ^ 
^c{t ^*)5aj,Pss^a /Tr(SA,Pss6A )■ This is the time-dependent 
density matrix conditioned on the detection of a red photon at 
Af at time t. The Liouvillian C is that of the free evolution of 
the total system. In the state p{t'), entanglement between the 
remote oscillators can occur 

With our assumptions about the cavity and oscillator pa- 
rameters, the effective phonon number rtM and linewidth 7 of 
the two oscillators will be approximately equal, such that we 
can drop indices on these quantities. For r ^ k~^,X~^, we 
find that IMI 



■) 

-AA 



(5) 



where the phase </> depends on path length differences. One 
could make this phase adjustable by introducing a phase shift 



in one of the arms of the interferometer, as suggested in Fig.[T] 
The interference effect in Eq. (|5]l can be understood classi- 
cally for large tim- It simply means that the red light phase 
difference between the two cavity outputs at time t is related 
to the blue light phase difference att + r because the mechan- 
ical oscillators have a well defined phase difference for times 
T < 7~^. This phase difference becomes time dependent if 
the mechanical frequency difference 5 is nonzero. From a 
quantum mechanical point of view, the detection of a red pho- 
ton creates a mechanical superposition with a definite phase 
6. The subsequent blue photon will be superposed between 
the upper and lower arms of the interferometer, with a relative 
phase determined by 6. This again determines the detection 
probability at Ab or Bb. A nonzero S means that the super- 
position switches back and forth from symmetric to antisym- 
metric as T increases, which is observable if 5 > 7. Note also 
that the ratio 5^ /(^AbiA ~'"3Bb|A ) can exceed its classical 
bound of 2/3 and come close to unity for small mean phonon 
numbers um- This means that, in the limit nutS/j 0, 
one can e.g. arrange the phases in such a way that when a red 
photon is detected at Ar, the next blue always arrives at Ab. 

We now discuss how entanglement between the two me- 
chanical oscillators can be detected. We wish to verify that 
following a red photon detection, the subsequent blue pho- 
ton is in a superposition state of originating from cavity 1 and 
cavity 2, i.e., that there is entanglement between the output 
modes 6b, 1 and 6b. 2- This entanglement must be smaller than 
or equal to the entanglement between the mechanical oscilla- 
tors. Again, we let t' = t + t. If p{t') is a separable state, it 
is straightforward to show that 



R{t) 



(S^,l(i')^b.l(i')^L(i')^b,2(i'))A. 



|(&Li(i')6b,2(t')) 



> 1 



(6) 
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This also holds for classical correlations between the two 
fields, which e.g. could originate from technical laser noise. 
Thus, we can use R{t) as an entanglement witness 1 18, 34 1. 

Expectation values with respect to the state p{t') can 
be measured through degrees of higher-order coherence, as 
Eq. (j4| suggests. We find that R{t) < R^{t), where the 
measurable upper bound for the entanglement witness is |34 | 



i?m(r) 



9aI\aM) ^ 9b!\aM) 



(2) 



(7) 



for a symmetric setup where the output flux is the same from 
both cavities. Thus, a measurement of Rm{T) < 1 is evidence 
of entanglement between the mechanical oscillators. 

We now show that the separability criterion ^ is violated 
as the oscillators are cooled close to the ground state. Inserting 
the expressions in Eq. (j5]l, we arrive at 



i?m(T) 



4nM [^M + {nu + l)e '^'^j 
(um + l)2e-27r cos2(,5t + 20) ' 



(8) 



againforr^K ^. In Fig.[2]we plot max(l — i?jjj(r), 0) 
as a function of phonon number um and delay time r in 
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units of 7~^, when assuming 6t + 2(f) to be an integer of tt. 
We observe that entanglement can be verified through viola- 
tion of the separability criterion for mean phonon numbers 
riM < 0.26. For um <^ 1, the entanglement can be detected 
for times r < ln((-\/2 — l)/2nM). It is also worth men- 




FIG. 2: (Color online) Density plot of max(l — R^{t),0) as a func- 
tion of mean phonon number um and time delay r . We set 5r + 20 = 
0. The separability criterion is violated below the white dashed 
curve, proving entanglement between the two mechanical oscilla- 
tors. The entanglement can be verified for mean phonon numbers 
Um < 0.26 and time delays r < In((\/2 — l)(nM + l)/2nM). 



tioning that before the first blue photon is emitted, the entan- 
glement between the mechanical oscillators only decays on 
the timescale of the mechanical phase decoherence, which is 
likely to be much larger than 7^^. This could be measurable 



by looking at the statistics of only the first blue photon follow- 
ing a red. 

To show that our proposed experiment is within reach 
of present day experiments, we take the membrane-in-the- 
middle geometry 1, 16,1 as an example. A set of realistic pa- 
rameters is wm/Ztt = 2 MHz, ujm/j ~ 2 • 10'^, K/27r = 1 
MHz, and |a|/27r = 10 kHz. Assuming an initial temper- 
ature of 20 mK, this gives an effective mean phonon number 
riM = 0.068, with riopt — 0.016. The output flux of red (blue) 
photons would be 41 (172) photons per second. The separa- 
bility criterion would be violated for times r < 0.47 millisec- 
onds. Note that this is over 400 times longer than the entan- 
glement lifetimes reported in a corresponding experiment with 
atomic ensembles by Chou et al. [24|. The long decoherence 
time is a result of high mechanical quality factors even when 
laser cooling to the quantum regime. 

In conclusion, we have proposed an experiment for entan- 
gling remote mechanical oscillators. This would be an impor- 
tant milestone in the endeavour to explore quantum effects in 
macroscopic systems. The entanglement is induced by opti- 
cal measurements and can be verified through second-order 
coherences of the optical field. Our proposal is relevant to 
present day experimental setups. We estimate entanglement 
storage lifetimes of milliseconds for the membrane-in-the- 
middle geometry, which could be of technological interest. 
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Supplementary Material to "Proposal for 
entangling remote micromechanical 
oscillators via optical measurements" 



I. INTRODUCTION 

This is a supplementary to the Letter "Proposal for entan- 
gling remote micromechanical oscillators via optical measure- 
ments". In Section [n] we briefly point out why our proposed 
experiment can also be relevant for systems that are not in 
the resolved sideband regime. Section |lll] provides some de- 
tails on how to calculate the correlation functions presented 
in the Letter In Section [rV] the challenge of filtering out the 
sidebands into separate spatial modes is addressed quantita- 
tively, whereas Section |V]presents a specific example of such 
a filter. Most importantly, in Section[Vl]we show how the de- 
grees of second-order coherence discussed in the Letter can be 
obtained in heterodyne photodetection without the need for a 



physical frequency filter. Finally, Section VII presents details 
on how our separability criterion can be related to observable 
correlation functions, and how it connects to a more familiar 
entanglement measure, namely concurrence. 



II. RESOLVED VS UNRESOLVED SIDEBAND REGIME 

In the Letter, we focus on the situation where the tempera- 
ture of the bath is large compared to the mechanical frequency, 
such that nth ^ 1- The oscillator then needs to be laser cooled 
by detuning the laser to A = —ujm- We consider the regime 
where the optical contribution to the mechanical damping 7opt 
is large compared to the intrinsic damping rate 7 to the ther- 
mal bath. In that limit, the effective phonon number is 



7 
7opt 



nth + nopt 



(9) 



To reach the quantum regime um < 1, one must then be 
in the resolved sideband regime ojm > k, such that nopt — 
(k/4wm)^ is small. 

However, one can also imagine a scenario where nopt < 1 
is not necessary to obtain entanglement. The oscillator would 
then need to be coupled to a very cold thermal bath, such that 
nth ^ 1- This is particularly relevant to optomechanical 
realizations with large mechanical frequencies. If one then 
operates in a parameter regime where the intrinsic damping 
rate 7 is large compared to the optically induced damping rate 
7opt (e g- through weak driving), the effective phonon number 
would be 

7opt 
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(10) 



which could still be small enough to achieve entanglement. 

III. DETAILS OF THE CALCULATION 

In this section, we present some details on how to calculate 
the correlation functions discussed in the Letter. 



A. Single cavity 

We first consider the case of a single cavity. The output field 
at the right-hand side of the cavity in a frame rotating at the 
drive frequency oj^ is b{t) e''^°*fe(t). Its Fourier transform 
is defined as 



We also define 



and the spectra 



die*"* 



dte"^*{b\t)b{0)) 



and 



(11) 



(12) 



(13) 



(14) 



These spectra can be found by solving the quantum Langevin 
equations in Fourier space and using the properties of the 
noise operators Cm and £^i. Let us define the bare cavity sus- 
ceptibility as 



1 



K/2-i{uj + A) ' 



(15) 



where A = ujj^ — ujc — g{x) is the detuning. We also define 
the effective mechanical susceptibility 



1 



7/2 - i{uj - (Dm) 



(16) 



For the spectrum flM, we have 6*^1 ^[w] — KFi|ap27r5(w) + 
kr I a \'^x~pi I xc [— wjp S'^ s [w] , where the mechanical oscillator 
spectral density is approximately 



IxmHP (nu + 1) + Ixh^^lP^-M 



(17) 



in the weak coupling Umit. 

The output filtered around the blue (red) sideband is defined 

by &b(r)(<) = e-'"°*6b(r)W with 



kit) = 



2tt 



and 



1 

2^ 



du e 



*F,[uj]b[uj] 



(18) 



(19) 



The filter functions ^b(i) ["^1 will naturally depend on the spe- 
cific details of the filter They should have the properties that 
Fh(i) [uj] ^ 1 for frequencies a; in a range of size > 7 around 
zbwM- On the other hand, |i^h(i)['^]| should be as small as 
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possible for a; ~ 0, such that the carrier photons at the drive 
frequency uju are rejected. The quantitative requirements for 
these filter functions are discussed in Section |IV] In the fol- 
lowing, we will assume that the filters are good enough to 
only pick out the sidebands, such that we can neglect the car- 
rier entirely. 

Since we linearize the equations of motion and only con- 
sider Gaussian input. Wick's theorem allows all higher-order 
correlation functions to be written in terms of correlation func- 
tions containing two operators. As an example, consider 



(2), ^ _ {bmbUt + T)b^{t + T)b,{t)) 



= 1 



{Umit)){bi{t)kit)) 

\m)k{t + T))\^ + \{bUt)blit + r))\ 



(20) 



{bim{t)){biit)b^it)) 



We will repeatedly make use of the fact that the effective me- 
chanical linewidth 7 is much smaller than the cavity linewidth 
K, as well as the integral 



duj 
2^ 



(21) 



valid for r > 0. We start with the red photon flux, which is 
given by 



ibim-it)) = 



cLu 
2^ 



(22) 



'«R|anXc[-WM]P('^M + 1) 
Similarly, the blue photon flux becomes 



2-K 



\F^,[-u:]\'S~,rM (23) 



= KR|a|^|xc[wM]|^ riM 
Moving on to the cross-correlators, we find 

duj 



{h\{t)bl{t + T))^ 



2-K 



2ti 



7 IXM[w]P(nM + 1) - XmM 



(24) 



We assume that the delay time r ^ A ^, k ^, in which case 
the last term does not contribute, giving 

(5t(t)fet(i + ^)) (25) 

The filter functions might give rise to the phase factor e^*'', 
but this does not contribute in Eq. (20 1. Finally, we have 



{blit)bl,it + T)) = 



27r 







(26) 



Inserting these results into Eq. pO| gives the expression for 



B. Two cavitites 



We now move on to the two-cavity setup in Fig. 1 of the 
Letter. The calculation of the correlation functions 5Ab|A (^)' 

5i3b|A (^)' performed in a similar way as above. How- 

ever, one needs to relate the fields at the photodetectors to the 
cavity output fields via the beam splitter transfer matrix: 



&A(i) 
bB{t) 



1 

1 

71 



b2{t)e'* + ibiit) 
hit) + ib2it)e"f' 

,(2) 



(27) 



The second-order coherences g^bi^ (''') g^^^j^^ (t) then 
follow from the correlation functions calculated in the previ- 
ous subsection, as well as the fact that the fields bi and 62 are 
uncorrelated in the absence of measurements (i.e. in the state 

Pss)- 



IV. REQUIREMENTS FOR THE FREQUENCY FILTER 

If one would like to detect the red and blue photons individ- 
ually on separate photomultipliers, one needs to separate the 
optical sidebands at the frequency cjd ± from the carrier 
at uju- This is a big technical challenge, partly because the 
mechanical frequencies can be in the kHz-MHz range, and, 
most importantly, because the carrier photons far outnumber 
the sideband photons in the weak coupling regime we discuss. 
Let us consider the photon flux of the light leaking out at the 
right hand side of the cavity. We denote /b (/r) the flux of blue 
(red) photons and fc the flux of photons at the laser frequency 

WD- 

The photon fluxes are 



/b 

/r 
fc 



'^R|aPlxc[wM]P?^M , 
«R|aPlxc[-WM]P("M + 1) 



(28) 



The cavity susceptibility xci'^] is defined in Eq. (15 1. Re- 
membering that |a| = gxzpf |a|, this means that the ratio be- 
tween blue photon flux and carrier photon flux is 



/b .f9x^p{\ 



(29) 



for a detuning A — —ujm- The ratio between the red and blue 
photon flux is 



A _ "-opt("M + 1) 

fh ("-opt + l)n-M 



(30) 



According to Eq. (29i, the relative flux of the sidebands 



to the carrier depends on the ratio between the bare optome- 
chanical coupling gxzpf and the cavity linewidth k. We will 
continue with the membrane-in-the-middle setup as an exam- 
ple. In present-day membrane-in-the-middle experiments, the 
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bare coupling constant gx.^-p{/2TT ~ 10 Hz f2\. If we op- 
timistically assume that this could be increased to 100 Hz, 
still use k/2tt — 1 MHz, and assume riM = 0.1, we find 
/b//c ^ 10^®. This indicates how challenging the filtering 
can be, at least in the membrane-in-the-middle setup. 



V. THE FREQUENCY FILTER - AN EXAMPLE 

In this section, we present an example of a specific filter that 
would separate the sidebands into different spatial modes. We 
do not claim that this is the best way to perform the sideband 
filtering. However, we wish to argue that, despite the con- 
siderable technical challenge, the task can be within reach of 
experiments. 

We imagine a filter consisting of four stages, as depicted in 
Fig. [3] The first stage is to combine the output field b{t) with 
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FIG. 3: (Color online) A simple suggestion for how to filter out the 
sidebands from the carrier. The signal is combined with a local os- 
cillator on a beam splitter (BS) with high reflectivity to null out the 
carrier as much as possible through destructive interference. It then 
proceeds to filter cavity 1, which mostly transmits carrier photons 
while reflecting the red and blue sidebands. The quarter wave plate 
(QWP) and the polarizing beam splitter (PBS) ensure that none of the 
reflected photons are lost. Filter cavity 2 transmits the blue sideband, 
whereas cavity 3 transmits the red sideband. 

a local oscillator (split off from the original laser) on a beam 
splitter with very high reflectivity. By adjusting the power and 
phase of the local oscillator, one can remove a significant part 
of the carrier at this stage. Ref. |3 1 reported an extinction ratio 
of 1/700 by this method. We will assume that the carrier flux 
is reduced by the factor R by this displacement. 

The next stage of our filter is the optical cavity labeled "Fil- 
ter cavity 1" in Fig. |3] The purpose of this cavity is to have 
the sidebands be promptly reflected, whereas a large part of 
the carrier is transmitted through the cavity. Filter cavity 1 
has a mode at the frequency lji with the linewidth fi. We 
require Fd ^ M ^ with Fd being the laser linewidth. 
Also, we let uji — ojd — Ai with Ai <C /i. Finally, we de- 
note the decay rate through the upper mirror and assume 



that it is /iu = /^(l/2 + 5i), i.e. approximately half the total 
cavity decay rate, with a deviation (5i <C 1. We assume that 
the initial signal b{t) is polarized. The polarizing beam split- 
ter (PBS) and the quarter-wave plate (QWP) then ensure that 
whatever is reflected from the cavity is also reflected by the 
PBS afterwards, such that we (in the ideal case) do not lose 
any of the photons that are reflected. We neglect complica- 
tions that might arise from a polarization mismatch between 
the incoming field and the internal cavity field. The vacuum 
noise at the PBS can be neglected, since it does not produce 
photodetection events. 

The purpose of filter cavity 2 is to let the blue sideband be 
transmitted, but reflect the red sideband plus what is left of the 
carrier Filter cavity 2 has a mode at the frequency uj2 ~ + 
wm and a linewidth A with the requirement 7 ^ A <C lom- 
Ideally, it should also be symmetric and have minimal loss, 
such that the left and right decay rates are equal, i.e. Al ~ Ap. , 
and Al = A(l/2 + 52) with 82 < 1. Filter cavity 3 has 
the same properties as cavity 2, except that it is centered on 
the frequency W3 = — wm, such that it transmits the red 
sideband. 

The three filter cavities results in the filter functions (see 



Eqs. (18l,(19l) 



PhM = — n \P^^\ 

A/2 — I — ujm) 



(31) 



FM = 



a/AlAt 



\52 + — ujm) 



A/2 — i{uj + ujm) a/2 — i{uj — ujm) 
where we defined the function 

fiSi + i{ui + Al) 



m/2-«(w + Al) 



(32) 



After the displacement by the local oscillator, the original out- 
put spectrum in the rotating frame is 



(rD/2)^ + c.2 

7 



-^■^c + /r,^2 , f ~ n2 /r 

(7/2) +[uj-ujm) 



(7/2)' + {oj + wm)' 



(33) 



when we assume the laser to have a Lorentzian lineshape with 
a width Fd. After the filtering, the blue output spectrum be- 
comes 



(34) 



By integrating over the (suppressed) carrier peak, we get the 
flux of remaining carrier photons in the blue output. 



(b) 



1 

2^ 



(35) 



AlAr / Fd 



■'M 



45f i?/c 



Here, the subscript c on the integration sign denotes that only 
the carrier peak is integrated over We have also assumed that 



8 



Fd, a <^ ujm- With these approximations, we reach the same 
expression for fc , the remaining flux of carrier photons in 
the red output. 

Let us use the same numerical example as in the Letter, 
where WM/27r = 2 MHz and j/2t: w 0.4 kHz. We as- 
sume that the filter cavities have widths /i/27r = X/2tt ^ 
10 kHz. Furthermore, we assume that the laser linewidth 
is Fd = 0.1 kHz and that AS'f < 10"^. If we also as- 
sume that the initial displacement reduces the carrier photon 



we arrive at /i'^''//< 



10- 



givmg 



flux hy R = 10" 

/b//i^^ ~ 10^. This shows that the carrier can in principle be 
removed from the sidebands in this idealized setup. It should 
nevertheless be mentioned that total suppression of reflection 
from a cavity is in practice quite difficult, and a minimal re- 
flection of about 10% is typical [2J. 



VI. HETERODYNE PHOTODETECTION 

In this section, we point out that the degrees of second-order 
coherence discussed in the Letter are accessible through het- 
erodyne photodetection, where the blue and red sidebands can 
easily be distinguished in the Fourier domain with a digital fil- 
ter function of our choice. This means that the measurement- 
induced entanglement can be detected also in cases where fil- 
tering of the sidebands into separate spatial modes proves to 
be too difficult. 



A. Single detector 

First, let us discuss the case of a single cavity. Assume that 
the output from the cavity is combined with a strong local os- 
cillator on a beam splitter with reflectivity r. The transmitted 
part of the local oscillator and the reflected part of the signal 
beam are then detected at a photomultiplier. We let the local 
oscillator frequency be wlo — — ^w, detuned from the 
original drive frequency by the intermediate frequency ujip. 
This mixes the carrier frequency down to wip, and the side- 
bands to ojip ± ujm- We will denote the photocurrent operator 
by i{t) and its Fourier transform by From this, we can 
define 



«r(t) 

ih{t) 



duj 

cLu 
2Tr 



'F,[uj]i[ujiF+uj] (36) 



e-"^^Fh[io]i[ujiP +uj] 



(37) 



where the filter functions Fi-[uj] and F}^[uj] can be chosen 
freely. As before, we should make sure that Fb(i.)[a;] ^ 1 
over the width of the sidebands at u — ±lum, and that they 
have minimal support at w ^ 0. More specifically, it is neces- 
sary that |Fb[0]| « (/b//c)^'/'\ where /b (A) is the flux of 
blue (carrier) photons, and similarly for the red filter function. 

Note that Zr(t) and ih{t) are not directly measurable quan- 
titites, but are complex functions that can be calculated from 
measurements of the photocurrent i{t). We can also express 



these "sideband currents" as 

/•OO 



/oo 
dtiF,(t-ti)e'"-*iz(ti) , j=r,b (38) 
'OO 

where the filter functions in the time domain are 

/oo T 
-^e— (39) 
-oo 

The filter functions Fb and F^. should be chosen such that they 
are sufficiently localized both in frequency and time, meaning 
that Fj [uj] falls off quickly compared to the mechanical fre- 
quency UJM, whereas Fj{t) falls off quickly compared to the 
mechanical lifetime I/7. A suitable choice would be Gaus- 
sian functions. 



i^b(r)H , 



(40) 



with the width satisfying 7 <C A ^ wm- With the num- 
bers used above (uju/^n = 2 MHz, A/27r = 10 kHz), 
the filter functions at the drive frequency become -Fb(r) [0] ~ 
exp(— 2 • 10"*), which would fully suppress the carrier. In the 
time domain, the filter functions are 



^b(r)W 



2tt 



(41) 



which falls off fast enough that one can resolve photons on 
timescales much smaller than I/7. 

We start by calculating the expectation value 

dh / dt2e-^"'^(*i-*^) (42) 

-00 J —00 

xF:{t^h)F,{t~t2){i{ti)i{t2)) . 

By the same approach as in Ref. ([T], one can show that the 
photocurrent autocorrelation function is 

(titMh)) (43) 
= iGeyrj{rj{: I{t,)I{h) :) + (: /(ii) - ^2)) • 

Here, G is the photodetector gain, e is the electron charge, and 
?7 is the dimensionless photodetection efficiency. The colons 
indicate time and normal ordering |1|. The operator I{t) de- 
scribes the photon flux at the photodetector at time t and is 
given by 

/(i) - (l-r)IAol' (44) 

in the limit where the (classical) local oscillator field transmit- 
ted by the beamsplitter, \/\ — r|/3io|e*^i°, is strong compared 
to the output field from the cavity. The operator hit) is the 
cavity output field operator in the rotating frame. The phase 
B\o is dependent on the path-length difference between the lo- 
cal oscillator and the signal beam, but is not important for our 



purpose. The integrand in Eq. (42 1 has several terms that os- 
cillate at the frequency wjp. In the limit wip ^ wm, these 
terms can be neglected and one arrives at 



2^ 



+ r7^{bl{t)k{t))] , (45) 
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where t = I — r and 



B. Multiple detectors 



(46) 



is the time domain version of Eq. ( 19 1. 

The last term in Eq. ( [45] l is proportional to the flux of red 
sideband photons leaking out of the cavity, whereas the first 
term is a shot noise term. In practice, the way to get rid of the 
shot noise term is to integrate the photocurrent spectrum S[ll!] 
around the frequency = wjp — after subtracting the 
shot noise floor. We denote this integrated sideband weight 
by Wr and the shot noise floor by So = {Ge)^r]t\l3io\^. We 
then find that 



-A/2 



/tJr-A/2 



So) = rr^So{h\{t%{t)) 
In exactly the same way, one can show that 

{i{t)i^{t))=So(^^+rr^(hl{t)K{t)) 
with the blue sideband operator given by 

kit) = / dh F^{t - ti)b{h) . 



(47) 



(48) 



(49) 



The term containing the blue photon flux can also be deter- 
mined from the photocurrent spectrum by 



W^= du (^M - So) = r7^So{blit)hit)) , 

(50) 

where we have defined = wif + i^M- 
Next, we define the correlation functions 



AjfcW = {^j{t)^k{t + T)) 



(51) 



where j, k — r, b. These functions can be determined from 
measurements of the photocurrent i{t). We will assume that 
the delay time t ^ A^^, such that Fj{T) ~ 0. In that case, 
one can show that 



^jk{r) 



rr^So{b]{t)bk{t + T)) , 



(52) 



rrjSoe'^'^^^{b]{t)blit + T)) 



These correlation functions are all we need to determine the 



(2) 

degrees of second-order coherence g^i - (t), with i,j 



r, b 



and T ^ A ^, whenever p^s is a Gaussian state. Specifically, 
we have 



(2) 



(r) 



1 



|A.,(r)p + |r.,(r)|^ 
W,W, 



(53) 



Notice that it is not necessary to determine the parameters G, 
r and 77. Furthermore, the second-order coherences are not 
susceptible to a frequency-dependent gain, as long as the gain 
is approximately constant over a frequency window of width 
7- 



In the two-cavity setup of the Letter, we are also interested 

(2) 

in observing cross-correlators such as g^^^j^ (t). In that case, 
one has to perform heterodyne detection on both the fields 
bA{t) and 6b (t), which are the cavity output fields after having 
been combined at the beam splitter (see Fig. 1 of the Letter and 
Eq. (p7]i of Section IIIi. We assume that this is carried out in 



the way discussed above, giving two photocurrents ZA(i) and 
«B(i), which have spectra ^Aiw] and 6*3 [w], respectively. 
The integrated red sideband weight of the spectrum 5'a [w] 



Wa, = / duj {Sa[uj]- So) (54) 

Jlj,--A/2 

= r7jSo{bi^{t)bAAt)) , 

and similarly for the blue sideband as well as the sidebands 
of the spectrum 5*8 [i^]- As above, we can define correlation 
functions from the filtered "sideband currents" i^, i^, if, and 



Aa^bJt) ^ {^fHt)^'^{t + r)) 

= rvSo{bi^it)bB,it + T)) , 
rA,B.(r) = {^fHt)hHt + r)) 

From these correlation functions, the degrees of second-order 
coherence can be found. As an example, we have 

Although we have for simplicity assumed that the local oscil- 
lator strength is the same for both the A and B fields, it should 
be noted that the second-order coherences are not susceptible 
to differences in these. 



(55) 



VII. ENTANGLEMENT WITNESS 

In this section, we discuss the entanglement witness R{t), 
defined in Eq. (6) of the Letter. In Section | VII A| we show 
how it can be expressed in terms of measurable higher-order 



coherences of the electromagnetic field. Section VII B relates 
our separability criterion to an entanglement measure on the 
state of the two mechanical oscillators. 



A. Relation to degrees of higher-order coherence 

We start by noting that 



(5T^,(t')6b,l(t')^L2(*')&b,2(t'))A. 



(57) 



ffS|ba|A,(^) (&LiW^b,i(t))(&L2(i)&b,2(t)) , 



10 



where we have defined the degree of third-order coherence 
9mi(^) = with 
Nijk = Iljn=i,j.k{bln{t)bjn{t)) ■ Also, observe that 



i(6L(t')&b,2(i'))A,r 



(58) 



> 



(2) 
^BblA, 



Using the relation 2(6^ JO^A.W) = E.=i,2(^L W^b,, W), 
we can derive an upper limit for R{t) (defined in Eq. 6 of the 
Letter), given by 



i?m(r) 



^i'b,2|b,l|A, 



(ffilV(^)-5Bl|A.(^) 



(59) 



second-order coherence. Specifically, we have 

ffS|b.i|A,M = 5S|A.M+5S|A,(T)+5S|b,i(0)-2 

(60) 



5S|A..(-)+dA,M 



(2) 



1 . 



The last equality is a result of the lack of correlations between 
bh,i{t) and 5b. 2 (i) in the state psa, i.e. in the absence of mea- 
surements. Assuming that the blue photon flux is the same 
from both cavities, i.e. (^t, i(i)^b.i(i)) — (^b 2(0^b,2(0)' we 
have 



Then, Eq. (7) of the Letter follows from Eqs. ( 59 1, ( 60 1, and 



The numerator is not directly measurable in the setup in Fig. 1 
of the Letter. However, one can imagine a modified setup 
where it is measurable, as shown here in Fig. [4] In that case, 
the blue and red sidebands are separated before combining 
the output from the two cavities. This allows a measurement 

°^ 5'Db|Cb|A,(^) = 5K2|b,i|A,('^)' ^" evaluation of 

Rmir). Note that for practical purposes, and to avoid detec- 
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FIG. 4: (Color online) Schematic view of a proposed experimental 
setup that allows a direct measurement of 3^'^2|b i|a (''')■ filter 
removes photons at the drive frequency and separates the sidebands 
into different spatial modes. Combining the red sidebands (dashed) 
from the two cavities on a beam splitter can create entanglement be- 
tween the mechanical oscillators. The blue sidebands (dotted) are 
used to verify the entanglement. 



tion loopholes 14J, it might be better to do the measurements 
in two steps. The second-order coherences can be measured 
without the detectors Cb and Db, when replacing the upper- 
and lowermost beam splitters with mirrors. The third-order 
coherence can be measured without Ab and Bb, when remov- 
ing the upper- and lowermost beam splitters. 

The measurement of the degree of third-order coherence 
(7^^2|b i|A however not necessary. The reason is that 

in the case of Gaussian states p^s, all higher-order correlation 
functions g*^"' (r) with n > 3 can be expressed by degrees of 



B. Relation to concurrence 

For mean phonon numbers riM ^ 1, the mechanical os- 
cillators will for the most part occupy the two lowest energy 
states, |0) and If we truncate to this subspace of or 1 
phonon, the state of the two oscillators is described by a 4 x 4 
density matrix. Assuming that coherences between different 
total number of phonons can be neglected, the density matrix 
is 



(62) 



Here, pij is the probability of having i phonons in oscillator 
1 and i phonons in oscillator 2, whereas q is the coherence 
between one phonon being in either oscillator 1 or oscillator 
2, given by q = {c\ci). Since the problem is now reduced to 
measuring entanglement between two-level systems, we can 
use the concurrence as an entanglement measure 15] . For the 
density matrix (62 1, the concurrence is max(C(r), 0), where 



(Poo 
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q* 


Pia 





I 








Pll 



(63) 



No entanglement gives a concurrence of 0, whereas maximal 
entanglement is characterized by a concurrence of 1 . 

In the weak coupling limit we have discussed, where 7 <C 
K, we can relate the blue output mode 6b, i to the mechanical 
oscillator annihilation operator Ci by 

kAt) ~ -iKRae-'^°'xc[uJM,^]c^{t) (64) 
+ vacuum noise . 

This means that we can relate the entries of the density ma- 
trix (62 1 to expectation values of the operators 6b. « = 1, 2, 
since the vacuum noise will not contribute. As discussed in 
the Letter, we are interested in expectation values in the state 
conditioned on a red photon detected at Aj at time t. In that 
non-stationary state, the entries of the density matrix become 
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dependent on t = t' — t. The elements needed to evaluate the 
concurrence C{t) are 



9(r) 



Pii(t) = 



(KR|a||xc[wM]|)^ 

(KR|a||xc[wM]|)'' 



and 



(65) 



Poo(t) = 1+Pii(r) 

_ (foLl(^O^b,l(tO)A, + (&^(^')6b,2(tO)A, 

('«R|a||xc['^M]|)^ 



(66) 



For a separable state, the concurrence C (t) = 0, such that 



Poo(t)pii(t) 



> 1 



(67) 



must be fulfilled. By using poo < 1, Eqs. (65 1 and (67 1 pro- 
duces the separability criterion presented in Eq. (6) of the Let- 
ter. We emphasize that even though the discussion here was 
restricted to the subspace of or 1 phonon and other assump- 
tions, the separability criterion (6) in the Letter is of general 
validity. 
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